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Abstract 

We present the analysis of Suzaku observations of the young open cluster 
Westerlund 2, which is filled with diffuse X-ray emission. We found that the emission 
consists of three thermal components or two thermal and one non-thermal compo- 
nents. The upper limit of the energy flux of the non-thermal component is smaller 
than that in the TeV band observed with H.E.S.S. This may indicate that active 
particle acceleration has stopped in this cluster, and that the accelerated electrons 
have already cooled. The gamma-ray emission observed with H.E.S.S. is likely to 
come from high-energy protons, which hardly cool in contrast with electrons. Metal 
abundances of the diffuse X-ray gas may indicate the explosion of a massive star in 
the past. 

Key words: stars: winds, outflows - - ISM: cosmic rays - - ISM: individual 
(RCW 49) — ISM: supernova remnants — Galaxy: open clusters and associations: 
individual (Westerlund 2) 



1. Introduction 

Cosmic-ray particles with energies of < 10 15 eV are thought to be accelerated in our 
Milky Way. Supernova remnants (SNRs) have long been considered to be the main sources 
of the galactic cosmic-rays. In fact, their synchrotron emission reveals the existence of high- 
energy electrons in SNRs (Koyama et al. 1995), and their gamma-ray emission may indicate 
that of high-energy protons (Aharonian et al. 2004). However, it is not clear to what extent 
the isolated SNRs could account for the cosmic-ray particles in our Milky Way. 

Young open clusters are one of other candidates of the cosmic-ray sources. The HEGRA 
stereoscopic system of air Cerenkov telescopes has detected a gamma-ray signal inside the 
core of the OB association Cygnus OB2 (Aharonian et al. 2002). The Milagro gamma-ray 
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observatory found gamma-ray emission around the young open cluster Berkeley 87 (Abdo et 
al. 2007). These observations suggest that particles are accelerated around open clusters. 

Several models have been proposed as the origin of the gamma-ray emission from young 
open clusters. The gamma-rays could be produced through collisions of protons accelerated at 
shocks in massive stellar winds in the cluster (Giovannelli et al. 1996). Some of the protons may 
penetrate the winds of the massive stars toward the stellar surfaces and produce gamma-rays 
through hadronic interactions in the dense parts of the winds (Torres et al. 2004). High- 
energy leptons may contribute to the gamma-ray emission through inverse Compton scattering 
(Bednarek 2007; Manolakou et al. 2007). In addition to the stellar winds, pulsars in a cluster 
might be responsible for the gamma-ray emission (Bednarek 2003). The gamma-rays could 
also be produced through the photodisintegration of highly boosted nuclei followed by daughter 
deexcitation (Anchordoqui et al. 2007). 

Westerlund 2 is one of the clusters from which gamma-ray emission has been detected 
with H.E.S.S. (Aharonian et al. 2007). It is ionizing the large H II region RCW 49 (NGC 3247). 
The gamma-ray emission is extended (~0.2°), and the whole cluster is buried in it. Observations 
in the X-ray band are crucial to find the origin of the diffuse gamma-ray emission from the 
cluster and the mechanism of particle acceleration. In particular, the strength of non-thermal 
X-ray emission can be used to discriminate between the hadronic and the leptonic origin of the 
gamma-rays (Bednarek 2007), and estimate the time-scale of particle acceleration (Manolakou 
et al. 2007). In the X-ray band, although the cluster has been observed with Einstein (Hertz 
& Grindlay 1984; Goldwurm et al. 1987) and ROSAT (Belloni & Mereghetti 1994), strict 
constraint on non-thermal emission has not been obtained. With Chandra, Townsley et al. 
(2005) studied diffuse X-ray emission < 5' from the center of Westerlund 2 and found that the 
spectrum could be represented by both thermal and non-thermal emission. 

In this paper, we report the results from Suzaku observations of Westerlund 2. Since 
Suzaku has a large collecting area and low background (Mitsuda et al. 2007), it is the best 
instrument for observations of dim and diffuse X-ray emission. This paper is organized as 
follows. The observation is presented in section 2. In section 3, we explain the details of 
the spectral analysis. Section 4 is devoted to discussion. Conclusions are summarized in 
section 5. We assume that the distance to the cluster is 5.4 kpc unless otherwise mentioned 
(Furukawa et al. 2009), although there is a debate about it (~ 2-8 kpc; Piatti et al. 1998; Dame 
2007; Tsujimoto et al. 2007; Rauw et al. 2007; Ascenso et al. 2007). At a distance of 5.4 kpc, 
1' corresponds to 1.6 pc. 

2. Observations 

Westerlund 2 was observed with Suzaku on 2008 August 9-11 and 2009 February 4-5. 
The exposure times for the two observations are 73.7 ks and 33.5 ks, respectively. The data 
of the two observations are dealt together in the following analysis. Suzaku has three working 
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Fig. 1. XIS raw image of Westerlund 2. The image is not corrected for vignetting and background. 
The calibration sources at the corner of the field were excluded. 

XIS CCDs (Koyama et al. 2007). Two of them (XIS 0, and 3) are front-illuminated (FI) and 
one (XIS 1) is back-illuminated (BI). The XIS was operated in the normal full-frame clocking 
mode. The edit mode was 3x3 and 5x5 and we combined the data of both modes for our 
analysis. We employed the following calibration files: XIS (20090203) and XRT (20080709). 
Events with ASCA grades of 0, 2, 3, 4, and 6 were retained. We excluded the data obtained 
at the South Atlantic Anomaly, during Earth occultation, and at the low elevation angles from 
Earth rim of < 5°. We further excluded the data with elevation angle from the rim of the 
shining Earth smaller than 20°. Figure. 1 is a raw XIS image of Westerlund 2. 

3. Spectral Analysis 

We analyze the spectrum of the diffuse emission around Westerlund 2. Since the angular 
resolution of the XIS is moderate (~ 2'), we need to estimate the amount of X-ray emission from 
point sources that contaminate the diffuse X-ray emission from the cluster (e.g. Hamaguchi et 
al. 2007; Hyodo et al. 2008; Ezoe et al. 2009). Moreover, since the diffuse emission is faint, the 
background spectrum must be constructed carefully. Errors on fitted spectral parameters are 
given at the 90% confidence level from now on. 

3.1. Contamination from Point Sources 

In order to estimate the contamination of the point sources, we used Chandra archive 
data of Westerlund 2 (ID 3501, 6410, and 6411). In order to extract the positions and spectra 
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Fig. 2. Simulated XIS image of Westerlund 2. Only photons from point sources are included. The 
exposure time is assumed to be 22 Ms in order to emphasize the X-ray structure. 

of the point sources, we used the ACIS Extract (AE) software package 1 (Broos et al. 2002). 
The procedures used in AE are described in Townsley et al. (2003) and Getman et al. (2005). 
AE detected 930 point sources from the data of the three observations (^ 10 _15 ergss _1 cm -2 in 
the 2-10 keV band). AE produces spectral files for individual point sources. For the 9 brightest 
sources with photon counts of > 300, we can treat their spectra separately. For fainter point 
sources, it is difficult to obtain their individual spectra because of their small photon counts. 
Thus, we combined their spectral files by MATHPHA included in FTOOLS. For that purpose, we 
used the data of 570 point sources, which are well-separated from their neighboring sources, 
and their background spectra can be obtained around them. We assume that the spectral 
shapes of the individual fainter sources are the same as that of the combined spectrum but 
their luminosities are different. 

Using the positional and spectral data, we simulate an XIS observation of the point 
sources with the XIS simulator MKPHLIST and XISSIM (Ishisaki et al. 2007). The produced 
event files include both positional and energy information of photons only from the point 
sources. Figure 2 is the simulated XIS image of Westerlund 2 that includes only photons from 
the point sources. 

1 The ACIS Extract software package and User's Guide are available online at 
http: / /www. astro. psu.edu/xray/docs/TARA/ae_users_guide. html. 
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3.2. Background Spectrum 

Since the TeV gamma-ray emission from Westerlund 2 covers the entire XIS field 
(Aharonian et al. 2007), we used the spectrum of a blank region around SWIFT J1010. 1-5747 
obtained with Suzaku (10 h 10 m 55?4, -57°51 m 14 s [J2000.0]) as the background of Westerlund 2. 
The exposure time is 38 ks. We call this region the blank region. This region is ~ 1.7° away from 
Westerlund 2 (10 h 23 m 59?0, -57 d 44 m 40 s ). The Galactic coordinates of Westerlund 2 and the 
blank region are (l,b) = (284.25°, -0.32°) and (282.87°, -1.38°), respectively. Westerlund 2 and 
the blank region are both affected by the Galactic Ridge X-ray Emission (GRXE; Revnivtsev 
et al. 2006). 

We subtract the non-X-ray instrumental background (NXB) from the spectrum of the 
blank region. The NXB is constructed from night Earth data, and is generated by the rou- 
tine XISNXBGEN included FTOOLS, considering the time-variation of the NXB (Tawa et al. 
2008). We calculate the effective area for the Suzaku XIS chips using XISSIMARFGEN, which 
provides the ancillary response file (ARF) through Monte Carlo simulations (Ishisaki et al. 
2007). After the NXB is subtracted and removing a few noticeable point sources brighter than 
~ 10~ 13 ergss~ 1 cm -2 in the 2-10 keV band, we fit the X-ray emission from the blank region 
with spectral components of the GRXE, the cosmic X-ray background (CXB), and the local hot 
bubble (LHB). We assume that the contribution of faint point sources (^3x 10~ 15 ergss~ 1 cm -2 
in the 2-10 keV band), which Chandra could have detected with an exposure time of ~ 100 ks, 
to the Suzaku spectrum of the blank region can be ignored, because Ebisawa et al. (2005) 
found that the contribution of point sources brighter than ~3x 10~ 15 ergss _1 cm -2 (2-10 keV) 
to the GRXE is only ~ 10 % (see also Revnivtsev et al. 2009). The spectrum of the GRXE 
can be represented by two thermal models (VAPEC), and two absorptions (PHABS) with a 
combination of PHABS * VAPEC + PHABS * VAPEC. The spectra of the CXB and the LHB 
are represented by an absorbed power-law model (PHABS * POWER) and a thermal model 
(APEC), respectively. For the CXB spectrum, we use the results of Kushino et al. (2002) for 
POWER and those of Dickey & Lockman (1990) and Kalberla et al. (2005) for PHABS. Since 
the region around Westerlund 2 is affected by the GRXE, we use the average of the LHB spectra 
at positions ±5 ° away from Westerlund 2 in the direction of galactic latitude (Snowden et al. 
1998). The temperature of the LHB is 0.1 keV. The parameters for the GRXE in the blank 
region are shown in Table 1. The best-fit values for the GRXE are derived by fixing all the 
parameters of the spectra of the CXB and the LHB at their best-fit values. For the groupings 
of metals, we refer to the results in Ebisawa et al. (2005). For the errors, we include not only 
statistical errors but also systematic errors made by a spatial variation of the intensity of the 
CXB. Kushino et al. (2002) estimated the spatial variation of the CXB is 6.49 % for the ASCA 
field. We adjusted the value to be consistent with the XIS field. Moreover, since we masked 
the regions around bright point sources (~ 20% of the XIS field), we consider the decrease of 
the effective area in the field. Thus, we adopt the systematic error of 16% for the CXB flux. 
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Table 1. Best-fit parameters for the GRXE 



Component 


Parameter 






Soft 


kT (keV) 




41 +0 ' 28 




Abundance except for Ne, Mg, Si 


(solar) 


0.044* 




Ne (solar) 




22+ - 07 




Mg (solar) 




ll+ 09 

u - ii -0.08 




Si (solar) 




25+ ' 35 

U - ZO -0.25 




N H (10 22 cm- 2 ) 




lfi+ 013 

u - iD -0.16 




Observed flux' (10 _8 ergs cm~ 2 s 


- 1 sr- 1 ) 


1 7 +al 

-0.0 




Intrinsic flux' (10~ 8 ergs cm~ 2 s~ 


1 sr- 1 ) 


Q+0.1 


Hard 


kT (keV) 

Abundance except for Fe (solar) 

Fe (solar) 

N H (10 22 cm- 2 ) 




0.17* 
1 5+ 1 - 1 
7fi+ - 48 

u - ' D -0.41 




Observed flux' (10~ 8 ergs cm -2 s 


- 1 sr" 1 ) 


a 0+0.9 

D - J -0.7 




Intrinsic flux'' (10 -8 ergs cm -2 s~ 


1 sr- 1 ) 


8.01J; 2 




X 2 /d.o.f. 




353/331 



* Fixed at the values of the GRXE in Table 8 of Ebisawa et al. (2005). 
t In the 0.7-10 keV band. 



Considering the distance between Westerlund 2 and the blank region on the sky, we need 
to take account of spatial variation of the GRXE, when we use the spectrum of the blank region 
as the background spectrum of Westerlund 2. In the 3-20 keV band, the flux of the GRXE in 
the blank region is 2.0 ±0.3 x 10 _11 ergscm _2 s _1 deg _2 and is consistent with the value obtained 
with RXTE (2.3 ±0.4 x 10~ n ergs cm" 2 s" 1 deg" 2 ; Revnivtsev et al. 2006). At the position of 
Westerlund 2, the X-ray flux of the GRXE is 1.8 ± 0.4 x 10" 11 ergs cm -2 s -1 deg -2 (Revnivtsev 
et al. 2006), which is ~ 60% of that of the CXB. Thus, we change the normalization of the 
GRXE spectrum of the blank region accordingly, considering the systematic error of the flux. 
Moreover, we include the systematic error of the CXB flux. Since the area of the masked regions 
in the Westerlund 2 field is different from that in the blank region (see the next subsection), 
we adopt an error of 18%. 

3. 3. Diffuse Emission from Westerlund 2 

Including the leaked X-rays from the point sources and the X-ray background emission 
(GRXE, CXB, and LHB) estimated above, we analyze the XIS spectrum of the diffuse emission 
from Westerlund 2. We refer to the diffuse X-ray emission around Westerlund 2 excluding those 
unwanted components as the gas component. We calculate the effective area of the Suzaku XIS 
chips using XISSIMARFGEN. All spectra were rebinned to give a minimum of 20 raw counts per 
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spectral bin to allow x 2 statistics to be applied. The NXB spectrum calculated with XISNXBGEN 
is subtracted from the spectrum. In order to minimize the contamination from the point sources, 
we mask the central region of the cluster (< 2' from the center). We also mask 88 bright point 
sources with circles of 0.5'-2'. In total, we mask ~ 40% of the XIS field. 

In the spectral fits, FI and BI XIS spectra were respectively summed. We fix the 
parameters for the leaked stellar emission and the X-ray background except for the nor- 
malizations of the GRXE and the CXB spectra, which we allow to vary according to the 
systematic errors estimated in section 3.2. We find that the spectrum of the gas compo- 
nent can be represented by three thermal models or two thermal models and one power- 
law model with absorptions. That is, PHABS * VAPEC + PHABS * (VAPEC + VAPEC) or 
PHABS * VAPEC + PHABS * (VAPEC + POWER). We call the former 3T model and the lat- 
ter 2TP model. In each model, metal abundances are the same among thermal models with 
different temperatures. The abundances of N and O are linked to that of C, and the abundance 
of Ni is linked to that of Fe. The abundances of He, Al, Ar, and Ca are fixed at the solar values. 
The parameters of these models are common between the two types of XISs (FI and BI), except 
for the normalization of the highest temperature VAPEC model (3T) or the power-law model 
(2TP), which is the main component of the spectra. 

The results of the fits are shown in Figure 3 and Table 2. The temperatures of the 
thermal components for 3T model are kTi, kT 2 , and kT 3 . Their intrinsic fluxes in the 0.7- 
10 keV band are j\, f 2 , and f 3 , respectively. We assume that the X-ray emission comes from 
the entire XIS field (17'. 8 x 17'. 8), and we make up the X-rays from the masked regions with 
those from the surrounding regions. For 2TP model, kT 3 and f 3 are replaced by the photon 
index (r) and the non-thermal flux (/nt), respectively. The absorption Nhi is the one for the 
lower temperature component (kTi), and Nh 2 is the one for the higher. The errors in Table 2 
include both the statistical errors and the systematic errors from the uncertainties about the 
GRXE and the CXB (section 3.2). The low metal abundances compared with the solar values 
are similar to those of the Carina Nebula, which also contains massive stars (Hamaguchi et al. 
2007). It is to be noted that the contributing ratio of flux of the gas component, the leaked 
stellar emission, and the X-ray background to the total spectrum excluding the NXB in the 
2-10 keV band is 43:14:43 (see Figure 4). 

Naze et al. (2008) monitored bright stars in Westerlund 2 and found that X-ray luminosi- 
ties of several of the brightest stars vary about 30%. However, the sum of the luminosities of the 
9 brightest stars in Westerlund 2 is only ~ 30% of the total stellar luminosity of Westerlund 2. 
Moreover, it is unlikely that the flux variation of the stars synchronizes. Considering the fact 
that only 14% of the diffuse X-ray emission around Westerlund 2 is attributed to the leaked 
stellar emission, we can ignore the flux variation of the stars in the spectral analysis. 
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Fig. 3. X-ray spectrum of the diffuse emission from Westerlund 2 (crosses). The NXB spectrum is 
subtracted. The result of the fit is shown by the lines, while the lower panel plots the residuals divided 
by the 1 a errors. 
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Table 2. Best-fit parameters for the gas component 



Parameter 3T 2TP 



N H1 (10 22 cm- 2 ) 


0.96^3? 




1 i+O.O 
1 - i -0.5 




kTi (keV) 


0.12Toi 




u - iZ -o.oi 




N H2 (10 22 cm- 2 ) 






1 4+ - 2 




kT 2 (keV) 






1 o +0 - 3 




kT 3 (keV) 


4 2+ 2 ' 1 








r 






9 9+0.4 
z - z -0.4 




C, N, (solar) 


u - z ^-0.20 




0-57tJ; 2 3 




Ne (solar) 


0.04igS 




n n7+ 017 

u - u '-0.07 




Mg (solar) 


0.46t°;f 7 




n 77 +0.98 
u -' '-0.42 




Si (solar) 


0.35l°; 2 ° 








S (solar) 


0.86l°;« 




0-78^; 2 31 




Fe, Ni (solar) 






o+ - 33 




fi (ergscm~ 2 s^ 1 ) 


8.7±?;° x 10- 


-12 


1.3±g; 3 x 10- 


11 


/2 (ergscm~ 2 s _1 ) 


2.2t°oi x 10- 


-12 


3.0±°i x 10- 


■12 


/3 (ergs cm -2 s _1 ) 


4.6±o;l x 10 ~ 


12 






/nt (crgscrrr 2 s _1 ) 






5Atii x 10- 


■12 


X 2 /d.oi. 


1719.56/1725 


1729.06/1725 



4. Discussion 

^.1. Upper limit of Non-thermal Flux 

We found that the spectrum of the gas component of Westerlund 2 can be represented 
by three thermal models (3T) or two thermal models and one power- law model (2TP). Since 
the power-law component in 2TP model can be replaced by a thermal component, the power- 
law component when the contributions of the GRXE and the CXB to the total spectrum 
are minimum gives the upper limit of non-thermal X-ray flux from the cluster. Assuming 
that the non-thermal emission comes from the entire XIS field, the upper limit of the flux 
is / NT < 2.6 x Kr^ergscm^s- 1 (0.7-2 keV) and 6.1 x 10~ 12 ergs cm" 2 s" 1 (0.7-10 keV). In 
this estimation, we compensate for the X-rays from the masked regions with those from the 
surrounding regions. 

Although the photon counts are small for detailed spectral analysis, we found that X-ray 
surface brightness of a region that is > 8' away from the cluster in the XIS field (the center 
of the circle is shifted from the cluster center according to the surface brightness contours) is 
comparable to the one estimated from the blank region. Thus, the emission from the outside 
of the XIS field does not much contribute to the total diffuse emission from Westerlund 2. 
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Fig. 5. Multi-wavelength measurements of Westerlund 2 (radio: Whiteoak & Uchida 1997; Einstein: 
Goldwurm et al. 1987; Chandra: Townsley et al. 2005; Suzaku: this study [0.7-2 keV]; EGRET; Hartman 
et al. 1999; H.E.S.S.: Aharonian et al. 2007). Fermi one year flux sensitivity is shown by the thin solid 
line 3 . A theoretical prediction of Manolakou et al. (2007) is shown by the dotted line (their model of 
t = 10 5 yr). 

We also analyzed the HXD data of Suzaku. Adopting a systematic uncertainty of 5% 
for the NXB, we found that the upper limit of the X-ray flux from Westerlund 2 in the 15- 
40 keV band is 5 x 1CT 12 ergs cm -2 s _1 , excluding the CXB and the GRXE. We assumed that 
the spectrum of the CXB is represented by a power-law model multiplied by a HIGHECUT 
model in XSPEC. We adopted the power-law index of 1.29, the power-law normalization of 
9.412 x 1(T 3 photons cm 2 s 1 FOV x keV 1 , and the e-folding energy of 40 keV, following a 
standard model 2 . The spectrum of the GRXE is extrapolated from the one in the soft X-ray 
band. The upper limit of the flux from Westerlund 2 can be compared with the extrapolated 
values of 3T and 2TP models (^4x 10~ 12 ergs cm -2 s _1 ). Considering that the HXD has a 
wider field (34' x 34') than the XIS, this result indicates that it is unlikely that strong hard 
X-ray emission extends beyond the XIS field. 

4-2. Comparison with Previous Observations 

In Figure 5, we present multi-wavelength measurements of Westerlund 2. The radio flux 
is considered as an upper limit, because it could only partially include the non-thermal radio 
flux produced by the energetic electrons (Whiteoak & Uchida 1997). The Einstein observations 
probably overestimate the diffuse flux, because the limited spatial resolution may not allow the 
subtraction of the point sources well enough (Goldwurm et al. 1987). The flux obtained with 
EGRET may also need to be regarded as an upper limit, because of its low spatial resolution 

(~ n 

2 http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/pin_cxb. html 

3 http: / /www-glast. slac.stanford.edu/software/IS/glast Jat_performance.htm 
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With Chandra, Townsley et al. (2005) found that the diffuse X-ray emission from 
Westerlund 2 can be represented by three thermal components (kT = 0.1, 0.8, and 3.1 keV), 
which is consistent with our 3T model (Table 2). Townsley et al. (2005) indicated that the 
hottest component can be replaced with a power-law component (r = 2.3). This photon index 
is also consistent with that in Table 2. Townsley et al. (2005) estimated that the flux of the 
power-law component is 1.4 x 10~ 12 ergscm" 2 s" 1 (0.7-10 keV). The difference from our best-fit 
value (Jnt = 5.4 x 10~ 12 ergscm~ 2 s -1 ; see Table 2) could be due to the difference of the regions 
investigated, and the backgrounds adopted. Townsley et al. (2005) focused on a narrower region 
(< 5' from the center of Westerlund 2), and they took a background in the same Chandra field. 
Thus, the flux should be regarded as an lower limit. On the other hand, since the power-law 
component can be replaced by a thermal component, the flux could also be regarded as an 
upper limit (Figure 5). 

With the high-quality of the Suzaku spectrum and the wide region we covered, we 
think that we have obtained the strictest upper limit of the diffuse non-thermal emission from 
Westerlund 2 so far. 

4-3. Comparison with Theoretical Models of Particle Acceleration 

Since Westerlund 2 has more than ten O stars and two Wolf-Rayet (WR) stars (Rauw 
et al. 2007), the strong winds from them may currently accelerate particles, which may be 
responsible for the TeV gamma-ray emission from the cluster. Considering the source extention, 
it is unlikely that a single star produces the gamma-ray emission (Aharonian et al. 2007). 
Bednarek (2007) and Manolakou et al. (2007) calculated broad band spectra of Westerlund 2; 
they chose parameters so that the spectra are consistent with previous gamma-ray and X- 
ray observations. They predicted that synchrotron emission from accelerated electrons could 
be observed in the X-ray band. Assuming that the distance to Westerlund 2 is 8 kpc and 
particles are injected uniformly during the lifetime of a WR star, Bednarek (2007) predicted 
that the X-ray flux would be 1.3 x 10" 11 ergscm -2 s -1 at ~ 1 keV. Assuming that the distance 
to the cluster is 2.8 kpc, Manolakou et al. (2007) predicted that the X-ray flux would be ~ 2- 
3 x 10~ n ergs cm -2 s" 1 at ~ 1 keV, depending on the injection time-scale. On the other hand, 
our new Suzaku results show that the upper limit of the non-thermal emission is /nt < 2.6 x 
10~ 12 ergs cm -2 s _1 (0.7-2 keV), which is much smaller than the above predictions (Figure 5). 
It would be difficult to adjust the parameters of the theoretical models to match both the new 
X-ray constraint and the previous gamma-ray observations, because the models must satisfy 
the relatively large ratio of gamma-ray to X-ray energy flux. 

Recently, Fukui et al. (2009) discovered jet and arc-like molecular feature around 
Westerlund 2. In particular, the latter suggests a past stellar explosion (or past stellar ex- 
plosions) in the cluster. Based on this observation, we consider particle acceleration associated 
with the explosion. 
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The total mass of the molecular gas around Westerlund 2 is ~ 2 x 10 5 M & and the size 
is ~ 30 pc (Furukawa et al. 2009). If the molecular cloud before the explosion of the star had 
almost the same mass and size, the average proton number density is ~ 100 cm" 3 . If the star 
explodes ~5x 10 5 yrs ago in the molecular cloud with the kinetic energy of 10 51 erg, the current 
velocity and radius of the shock wave is ~ 16kms _1 and ~ 40 kpc, respectively (see equation 2 
in Yamazaki et al. 2006. The radius can be obtained by the integration of the velocity). The 
velocity is comparable to the internal velocity of the molecular gas, and the size is comparable 
to that of the arc (~ 30 pc) observed around Westerlund 2 (Fukui et al. 2009). Although the 
shock may still be expanding, it is unlikely that particles are accelerated at the shock at present 
because of the small velocity. However, particles might be accelerated in the past when the 
velocity of the shock was large. If adiabatic energy loss did not much affect the particle energy, 
the accelerated protons should not have lost their energy because of their long cooling time. 
Thus, the protons accelerated in the past may be producing the TeV gamma-rays through 
pp-interactions in the surrounding molecular gas. On the other hand, accelerated electrons 
should have lost their energy through synchrotron emission. Therefore, the current ratio of 
gamma-ray to X-ray energy flux should be large, because the electrons are the source of the 
X-ray flux through synchrotron emission (Yamazaki et al. 2006). This is consistent with the 
observations of Westerlund 2; the observed ratio of the 1-10 TeV energy flux to the 2-10 keV 
energy flux is i?TeV/x > 2.7, which actually indicates that the stellar explosion did not happen 
recently (Yamazaki et al. 2006). Such stellar explosions in dense clouds could be the sources 
of unidentified TeV sources and the so-called PAMELA anomaly (Bamba et al. 2009; Ioka & 
Meszaros 2009; Fujita et al. 2009). 

The mass of the star responsible for the explosion could be extremely large, because stars 
in an open cluster form almost simultaneously and massive stars with ~ 80 M & still survive 
in Westerlund 2 (Rauw et al. 2005). In general, such massive stars produce ejecta with a 
large a-element to iron abundance ratio (Kobayashi et al. 2006). This trend is consistent with 
the metal abundances shown in Table 2, although the contribution from stellar winds from 
surviving stars must be considered for detailed analysis. If the thermal gas in 3T model is 
uniformly distributed within 8' from the cluster center, the total gas mass is 890±^|8 M . From 
the metal abundances in Table 2, the masses of Si and S contained in the gas are estimated to 
be ~ 0.1-0.6 M Q and ~O.1-O.7M , respectively. Since a star with a mass of >4OM produces 
~ 0.3 M of Si and ~ 0.1 M of S (Kobayashi et al. 2006), explosions of a few massive stars are 
enough to produce observed Si and S, even if metals that originate from sources other than the 
stellar explosions are not considered. Since some of the stars with masses of ^ 30 M Q explode 
as hypernovae (Galama et al. 1998; Iwamoto et al. 1998), the one exploded in Westerlund 2 
might be such a hypernova and might trigger a gamma-ray burst. However, because of the 
large errors for the metal abundances in Table 2 and the lack of the knowledge about the mass 
of the progenitor star, we cannot firmly determine whether the progenitor star is a supernova 
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or a hypernova from the metal abundances (Kobayashi et al. 2006). 

It is to be noted that pure hadronic models could not account for the GeV gamma-ray 
radiation detected by EGRET at the position of Westerlund 2 (Hartman et al. 1999; Bednarek 
2007). In the near future, Fermi will reveal whether the GeV gamma-ray source is related to 
the TeV gamma-ray source in Westerlund 2 and by what mechanism the gamma-ray radiation 
is produced (Figure 5). 

5. Conclusion 

We observed the young open cluster Westerlund 2 with Suzaku X-ray satellite. We 
found that diffuse X-ray emission extends to ~ 8' from the cluster center. We analyze the 
spectrum considering the contamination from point sources using Chandra data. We found 
that the diffuse emission consists of three thermal components (kT ~ 0.1, 0.9, and 4 keV) or 
two thermal components (kT ~ 0.1 and 1.0 keV) and one non-thermal component (r ~ 2.2). 
The upper limit of the non-thermal energy flux is smaller than the TeV gamma-ray energy flux 
observed with H.E.S.S. The abundances of a-elements are relatively high in comparison with 
that of iron. 

The relatively high X-ray to gamma-ray energy flux ratio suggests that the gamma- 
ray emission is attributed to protons that were accelerated in the past, because the cooling 
time of the high-energy protons is much longer than that of high-energy electrons, and the 
electrons that are responsible for synchrotron X-ray emission may have cooled. Considering 
the structures in molecular gas surrounding the cluster, protons would have been accelerated at 
a shock created through a stellar explosion that happened ~ 10 5 -10 6 yrs ago. Since extremely 
massive stars (~ 80 M Q ) still survive in this cluster, the star responsible for the explosion would 
also have been very massive. The high abundances of a-elements compared with iron may 
support this idea. Thus, the progenitor star could have exploded as a hypernova rather than 
as a normal supernova. 

We wish to thank M. Tsujimoto for providing us the Chandra data. We also thank 
K. Hamaguchi, and Y. Hyodo for useful discussion. This work was supported in part by a Grant- 
in-aid from the Ministry of Education, Culture, Sports, Science, and Technology (MEXT) of 
Japan, No. 20540269 (Y. F.). 
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